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THE  RAYLEIGH-TAYLOR  AND  KELVIN-HELMHOLTZ  INSTABILITIES 
IN  TARGETS  ACCELERATED  BY  LASER  ABLATION 

I.  INTRODUCTION 

The  Rayleigh-Taylor  (R-TJ^instability  arises  in  the  acceleration  of  a 

fluid  by  one  of  lower  density.  This  instability  may  represent  an  obstacle 

to  inertial  confinement  fusion  by  destroying  the  spherical  symmetry  of 

high  aspect  ratio  imploding  shells.  The  R-T  instability  causes 

corregations  on  the  ablation  layer  between  the  cold,  high  density  shell 

and  the  hot,  low  density  ablated  material.  These  corregations  eventually 

2 

grow  nonlinear ly  to  form  the  typical  bubble-and-spike  structure.  In  addition 
the  shear  flow  that  develops  as  the  bubble-and-spike  structure  evolves  can 
lead  to  the  Kelvin-Helmholtz  (K-H) 3,4  instability. 

Although  much  research  has  been  applied  to  the  study  of  the  R-T  instability 
in  ablatively  accelerated  targets,  a  complete  analytical  theory  that  includes 
important  physical  phenomena  such  as  finite  density  gradients,  ablation 
effects,  mass  transfer,  convection  and  finite  thermal  conduction  is  very 
difficult  to  perform  and  still  does  not  exist.  The  nonlinear  growth, 
turbulence  and  saturation  processes  are  only  now  being  investigated  with 
computer  codes,  albeit  with  somewhat  conflicting  results.5'6,7 

II.  MODEL 

We  report  here  on  our  investigation  of  the  R-T,  K-H  instabilities  in 

8  9 

laser  ablatively  accelerated  flat  targets  using  the  FAST2D  '  laser-shell 
simulation  code.  This  is  a  fully  two-dimensional  Cartesian  code  with  a 
sliding  rezone  Eulerian  grid  with  variable  grid  spacing.  The  grid  spacing 
is  0.25  urn  for  10  zones  on  either  side  of  the  ablation  layer  and  increases 
uniformly  to  a  2.  um  spacing  for  most  of  the  rest  of  the  grid.  The  finely 
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zoned  region  near  the  ablation  layer  is  required  in  order  to  accurately  re¬ 
solve  the  steep  density  gradient.  The  zones  in  the  underdense  plasma  beyond 
the  critical  surface  and  in  the  low  density  rear  portion  of  the  foil  are 
stretched.  The  system  has  40  zones  transverse  to  the  laser  (y-direction)  and 
120  zones  parallel  (x-direction) .  The  system  is  periodic  in  the  transverse 
direction. 

FAST2D  solves  the  ideal  hydrodynamic  equations  using  the  flux-corrected 
transport  (FCT) 10  algorithms  with  a  two-dimensional  classical  (T^^2)  plasma 
thermal  conduction  routine.  FAST2D  accurately  models  rotational  flow  pheno¬ 
mena  -  a  necessary  criterion  for  any  numerical  code  that  intends  to  model 
the  R-T  instability  into  the  nonlinear  regime  with  its  accompanying  shear 
flow.  It  is  important  to  realize  that  rotational  flows  occur  with  the  R-T 
instability  irrespective  of  the  appearance  of  the  K-H  instability.11 
The  initial  density,  pressure  and  temperature  profiles  for  an  20  ym 

thick,  plastic  (CH)  foil  irradiated  with  an  absorbed  laser  intensity  of 
13  2 

1.0  x  10  Won  are  generated  from  a  one-dimensional,  analytic,  quasi- 

12 

static  equilibrium  model.  The  solutions  have  been  shown  to  have  provided 

an  adequate  steady-state  with  a  one-dimesional  laser-foil  computer  code.13 

A  representative  initial  profile  is  illustrated  in  Fig.  la.  The  initial 

conditions  are  assumed  uniform  in  the  transverse  direction  of  the  foil  . 

A  two-dimensional  perspective  plot  of  the  mass  density  is  shown  in  Fig.  lb. 

The  initial  perturbation  is  obtained  by  perturbing  the  density 

profile  at  its  peak.  The  first  ten  modes  are  excited  with  equal  amplitudes 

and  random  phases  corresponding  to  a  total  initial  density  perturbation  of 

8%  ( (P  -  P  )/<P>  =  0.08).  The  longest  wavelength  excited  is  100  ym. 
max  mm 
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LOG  IP. 


Fig.  la  —  A  semi-long  plot  of  the  initial  pressure,  density  and 
temperature  profiles  of  a  thin,  plastic  foil 


FA5T2D  LASER  SHELL  0.00*10  0 


Fig.  lb  —  A  two-dimensional  perspective  plot  of  the  initial  mass 
density.  The  laser  is  coming  in  from  the  front. 
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III.  NUMERICAL  RESULTS 


The  amplitudes  of  the  first  five  modes  are  shown  in  Fig.  2.  is  the 
Fourier  transform  of  the  summed  mass  from  the  rear  of  the  foil  to  the 
ablation  edge  for  each  y-coordinate.  Note  that  the  fifth  mode  [X^  =  20  um 
(the  foil  thickness) ]  does  not  undergo  any  serious  growth  until  after  7 
nanoseconds.  At  this  time  the  second  and  third  modes  undergo  a  slight  but 
noticeable  transition  in  their  growth  rate.  It  is  at  about  this  time  that  the 
bubble-and-spike  formation  becomes  quite  large  and  a  strong  shear  flow 


develops  which  accelerates  the  growth  of  the  spike.  The  growth  saturates 

at  about  9  nanoseconds  and  the  foil  fragments  at  10.6  nanoseconds. 

The  growth  rates  for  the  first  5  modes  are  compared  with  the  classical 

value,  y  =  (kg)  ,  on  Fig.  3.  These  growth  rates  are  well  below  the  classical 
c 

14 

values,  in  agreement  with  other  time  dependent  calculations.  We  also 
indicate  growth  rates  obtained  from  a  modified  solution  of  Hsieh's15  inves¬ 
tigation  of  the  R-T  instability  in  boiling  films.  This  solution  includes 
mass  and  thermal  transfer  across  the  fluid  layers.  The  growth  rate  can  be 
expressed  as 


,  .  a  .  2  a 

YH-  [<->  ♦  M  --  , 

where  p is  the  density  of  the  "heavy"  fluid  and 

G  ,1  1  , 

a  =  -  (r-  +  r~)  . 

L  hi  h2 

hl^h2^  t*le  depth  of  the  "light"  ("heavy")  fluid,  G  is  the  thermal  flux 
across  the  interface  and  L  is  the  latent  heat.  L  is  determined  by  calculating 
the  change  in  entropy  across  the  ablation  layer. 

Fig.  4  shows  a  time  sequence  of  density  contours  when  the  system  is  far  into 
the  nonlinear  regime.  The  corresponding  mass  density  perspective  plots  are  shown 
in  Fig.  5.  The  lines  are  contours  of  constant  density  in  10%  increments  of  the 
maximum  density  (0^^  *  0.71  gm/cm3)  counting  from  the  outside  inward.  By  8  ns 


Fig.  4  —  Density  contours  at  different  times  showing  the  development  of  the 
bubble-and-spike  formation  and  the  subsequent  tip-spreading.  The  rightmost 
contour  is  the  critical  surface. 


the  bubble-and-spike  structure  is  quite  prominent.  The  spikes  grow  rapidly  in 
amplitude  and  the  classical  tip  spreading  is  well  developed  at  the  shoulders  of 
the  spikes  by  9.2  ns.  As  the  spikes  accelerate,  lateral  flow  continues  to  widen 
the  tips  until  the  foil  fragments  at  10.6  ns. 

The  tipsreading  is  due  to  the  nonlinear  phase  of  the  Kelvin-Helmholtz  insta¬ 
bility  in  the  shear  layer  at  the  interface  between  the  falling  spike  and  the  rising 
bubble.  That  a  very  strong  shear  layer  develops  between  the  bubble  and  the  spike 
is  evident  in  the  vector  fluid  velocity  field  (Fig.  6)  spanning  the  region  between 
the  arrows  on  the  density  contour  plots  of  Fig.  4.  The  vector  lengths  are  scaled 
to  the  maximum  velocity  (VELMAX)  as  indicated  on  each  figure. 

The  lateral  flow  that  causes  the  tip  widening  is  clearly  apparent  in 
the  velocity  plots.  Also  note  the  strong  circulatory  flow  (vortex  flow) 
that  develops  on  the  spike  side  of  the  shear  interface.  Although  no  spiral 
roll-up  appears  in  the  wake  of  the  bulge,  the  existence  of  vortex  motion 

is  indicated  by  the  non-curl  free  nature  of  the  fluid  flow  at  the  spike 

shoulders  and  the  flattening  of  the  rear  of  the  bulge.3 

The  time  evolution  of  the  bubble-and-spike  formation  is  shown  in 
Fig.  7.  Up  to  about  10  ns  the  growth  of  the  spike  outpaces  that  of  the 

bubble.  Initially  the  spike  experiences  less  drag  than  the  bubble  since 

the  "heavy"  spike  is  replacing  fluid  with  a  smaller  density  whereas  the 
"light"  bubble  is  attempting  to  rise  through  a  denser  medium.  It  is  not 
until  the  area  of  the  bubble  becomes  quite  large,  thereby  increasing  its 
bouyancy,  that  the  bubble  velocity  approaches  the  spike  velocity. 

The  amplitude  of  the  bubble-and-spike  formation  (AX^)  is  plotted  as 

a  function  of  time  in  Fig.  8.  Here  we  have  also  plotted  the  distance  from 

the  critical  surface  to  the  spike  tip  (AX  )  and  the  distance  from  the  rear 

cs 

of  the  foil  to  the  bubble  (AX  ,  )  -  a  measure  of  how  far  the  bubble  has 

rb 

"risen".  The  spike  "falls"  nearly  exponentially  up  to  saturation  at  about 
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Fig.  6  —  Vector  fluid  velocity  plots  for  the  region  spanned  by  the  vertical  arrows 
of  Fig.  4.  In  each  figure,  the  vector  length  is  scaled  to  the  maximum  velocity 
(VELMAX).  Note  the  vortices  developing  at  the  bubble-and-spike  interface  at  9.2 
and  10  ns. 


Fig.  7  —  Bubble-and-spike  contours  at  various  times.  Note  the  large  increase 
in  the  area  of  the  bubble  between  9.33  and  10.01  ns. 
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Fig.  8  —  A  plot  of  the  spike  amplitude  (AX^),  the  distance  between  the  critical  surface 
and  spike  tip  (A  X^)  and  the  distance  between  the  rear  of  the  foil  and  the  bubble  (AX^) 
as  a  function  of  time 
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10  ns.  Note  the  reduction  in  the  spike  growth  rate  after  8  ns.  It  is  at 
about  this  time  that  the  tip-widening  appears.  The  bubble  "rise"  is  much 
slower  and  approaches  zero  between  9  and  10  ns.  During  this  period  the 
area  of  the  bubble  is  increasing  radically  (see  Fig.  7)  until  the  bubble 
bouyancy  is  large  enough,  coupled  with  the  reduced  mass  of  the  shell  due 
to  ablation,  that  the  bubble  bursts  through  the  rear  of  the  foil. 

The  temperature  contours,  Fig.  9,  indicate  a  strong  temperature  gradient 
at  the  spike  tip,  corresponding  to  a  large  mass-ablation  rate  there;  albeit, 
not  large  enough  to  cause  a  recession  of  the  spike  growth.7  Note  that  the 
temperature  of  the  bubble  is  quite  cold  and  uniform,  thus  mass  ablation  must 
be  removed  from  consideration  as  a  possible  stabilizing  mechanism  for 
the  K-H  instability  that  appears  at  the  spike-bubble  interface. 

The  K-H  instability  is,  in  fact,  the  stabilizing  mechanism  here.  The 
spike  tip  widening  is  due  to  the  growth  of  the  K-H  instability  in  the  shear 
layer  at  the  bubble-and- spike  interface.  This  tip  widening  increases  the 
frontal  area  of  the  spike,  thereby  increasing  its  drag  and  reducing  its 
rate  of  "fall".  With  a  thicker  shell  the  spike  tips  may  widen  enough  to 
seal  off  the  bubble  and  form  a  uniform  ablation  front  thereby  saturating 
the  spike  growth  well  before  shell  fragmentation.  This  aspect  of  the 
problem  is  now  being  investigated  and  is  the  subject  of  a  future  report. 

Since  the  growth  rates  of  the  instability  are  well  below  the  classical 
values/  at  10.  ns  the  rear  of  the  foil  is  still  being  accelerated  uniformly, 
as  shown  by  the  pressure  contour  in  Fig.  10.  At  this  time  the  foil  velocity 
is  160  km/s  -  within  the  range  required  for  inertial  fusion.  The  foil  has 
been  pushed  a  distance  of  800  urn  -  forty  times  its  original  thickness. 
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Fig.  10  —  Pressure  contours  at  10  ns.  Note  that  the  three  leftmost  contours  are 
parallel  indicating  a  uniform  acceleration  of  the  rear  of  the  foil. 


IV.  SUMMARY 


In  summary,  using  the  FAST2D  laser-shell  simulation  model  we  have 
accelerated  a  20  ym  thick  plastic  foil  up  to  160  km/s.  The  foil  maintained 
its  integrity  up  to  10.  ns  giving  an  aspect  ratio  of  40.  We  are  able  to 
follow  the  Rayleigh-Taylor  bubble-and-spike  development  far  into  the  non¬ 
linear  regime  and  beyond  the  point  of  foil  fragmentation.  Strong  shear 
flow  develops  at  the  interface  between  the  bubble-and-spike  which  leads 
to  the  development  of  the  Kelvin-Helmholtz  instability.  The  K-H  insta¬ 
bility  causes  the  tips  of  the  spikes  to  widen  and  as  a  result  reduce 
their  rate  of  "fall".  This  tip  widening  is  a  possible  stabilization 
mechanism  for  the  R-T  instability  for  a  thicker  foil  in  that  the  spike 
tips  may  widen  enough  to  seal  off  the  bubble  and  form  a  uniform  ablation 
front. 

The  nonlinear  aspects  of  these  results  differ  from  those  of  Ref.  7.  We 
believe  there  may  be  several  reasons  for  this  difference.  The  laser  intensities 
used  by  Ref.  7  are  two  orders  of  magnitude  larger  than  the  intensities  used 
here.  This  would  result  in  a  much  larger  mass  ablation  rate  and  could 
possibly  explain  the  spike  erosion  that  they  cite.  However,  since  the 
distance  from  ablation  to  critical  scales  as  the  laser  intensity,  their 
ablation  layer  should  be  on  the  order  of  3-5  mm  away  from  the  critical 
surface.  Coupled  with  this  large  separation  between  the  ablation  and  cri¬ 
tical  surfaces  is  a  temperature  gradient  scalelength  at  the  ablation 

16 

surface  on  the  order  of  6-12  ym  for  a  3  ym  carbon  slab  and  a  laser 
15  2 

irradiance  of  10  W/cm  .  It  is  difficult  to  understand  how  the  mass 
ablation  rate  could  be  sensitive  to  a  spike  liss  than  a  micron  in  length. 
Another  reason  may  stem  from  the  difference  in  the  numerical  algorithms. 17 
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